Electrically driven convection in a thin annular film 
undergoing circular Couette flow 
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We investigate the linear stability of a thin, suspended, annular film of conducting fluid with 
a voltage difference applied between its inner and outer edges. For a sufhciently large voltage, 
such a film is unstable to radially-driven electroconvection due to charges which develop on its 
free surfaces. The film can also be subjected to a Couette shear by rotating its inner edge. This 
combination is experimentally realized using films of smectic A liquid crystals. In the absence of 
shear, the convective flow consists of a stationary, azimuthally one-dimensional pattern of symmet- 
ric, counter-rotating vortex pairs. When Couette flow is applied, an azimuthally traveling pattern 
results. When viewed in a co-rotating frame, the traveling pattern consists of pairs of asymmetric 
vortices. We calculate the neutral stability boundary for arbitrary radius ratio a and Reynolds 
number TZe of the shear flow, and obtain the critical control parameter TZc {a, TZe) and the critical 
azimuthal mode number TOc {a, TZe). The Couette flow suppresses the onset of electroconvection, so 
that TZc {a, TZe) > TZc (a, 0). The calculated suppression is compared with experiments performed 
at Q = 0.56 and < TZe < 0.22. 
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I. INTRODUCTION 

We examine the effect of a circular Couette shear on 
radially driven convection in a two-dimensional annu- 
lar fluid. Surprizingly, this remarkably ideal situation 
is experimentally realizcable in the electroconvection of 
thin, freely suspended liquid crystal films. 0] This sys- 
tem, which can be accurately described by electrohydro- 
dynamic theory presents a unique opportunity to 

quantitatively study convection with a novel combination 
of forcing and shear. In this paper, our treatment is both 
theoretical and experimental; we present a complete lin- 
ear stability analysis and use it to make the first quantita- 
tive comparisons with experiment. [0 This rather simple 
instablity has a rich phenomenolgy in its nonlinear be- 
haviour, and is a promising new testing-ground for the- 
ories of bifurcations in nonlinear patterns. Q| Since it 
involves radial driving forces, this system may also be in- 
teresting as an experimental analog for some geophysical 
instabihties. [^,|| 

The linear instability and subsequent nonlinear evolu- 
tion of flows depend strongly on the symmetry and struc- 
ture of the unstable base flow. One way of systematically 
studying this dependence is to superpose simple flows or 
rotations on well-understood instabilities. For example, 
the standard case of buoyancy-driven Rayleigh-Benard 
convection (RBC) has been studied in the presence of 
rotation |^,p|,ppiC|] and shearing due to an open through- 
flow [pl] , pL2|[ 7 as well as in the geophysically interesting 
cases of radial gravitation with rotation [^|J|]. The phe- 
nomenology of two-dimensional (2D) electroconvection 
in a rectangular geometry was the subject of our pre- 
vious experimental |l^,|l^,|l^,|l6|,0 and theoretical 
work. These studies made precise the degree of analogy 



with RBC; the electrical nature of the forcing introduces 
some crucial differences in detail, even at the level of 
linear stability. Here we superpose this basic instabil- 
ity on a shear in the form of a circular Couette flow. 
The annular geometry and the Couette shear naturally 
bring to mind another standard instability, that of Cou- 
ette flow between concentric cylinders, which is unstable 
to three-dimensional (3D) Taylor Vortex flow (TVF). 
TVF has also been examined for the cases of an imposed 
axial throu gh-fl ow |18| and superposed on radial electri- 
cal driving |]l9| . Although our analysis involves a similar 
cylindrical geometry, the crucial difference in our case is 
that the flow is truly 2D, which prevents the 3D Taylor 
instability. ||2^,^ We will show, in fact, that the Couette 
shear is always stabilizing in our system. 

In summary, the unique fluid-mechanical features of 
our system are that it achieves the superposition of a 
circular Couette shear with a radially (transverse to the 
shear) driven convective instability in an entirely 2D flow. 
This combination cannot be accurately realized in any 
other experimentally accessible system. We will show 
that this combination leads to some interesting new phe- 
nomena. 

Electroconvection results when an applied electric 
field acts on space charges within a fluid, resulting in 
body forces which drive flows. Electroconvection in 3D 
isotropic dielectric fluids has been the subject of many 
previous studies. |^2| In general, the agreement between 
experiment and theory has been modest, owing to diffi- 
culties in modeling the generation of charge at the elec- 
trodes. 

Another often-studied form of electroconvection occurs 
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FIG. 1. Schematic of an annular film and cylindrical coor- 
dinate system. The film has a thickness s ^ d in the z di- 
rection (out of the page), and is treated as a two-dimensional 
sheet. 

in certain anisotropic 3D fluids, i.e. liquid crystals, in 
which bulk charge separation occurs due to a mechanism 
which is unique to anisotropic fluid conductors, Since 
this mechanism does not involve electrode effects, signifi- 
cant agreement between theory and experiment has been 
achieved despite the complexity of the anisotropic fluid. 

In thin, suspended fluid films, another mechanism of 
charge separation can take place, leading to 2D electro- 
convection. If a voltage is applied across a conducting 
fluid film, the free surfaces acquire a potentially unsta- 
ble distribution of surface charge. This charge must 
be present simply to satisfy the standard matching con- 
ditions on electric fields across an interface. We do not 
consider charge generation by electrode injection or other 
bulk effects. A non-zero threshold voltage for the onset of 
electroconvection results from the competition between 
the electrical force on the surface charge configuration 
and viscous and electrical dissipation. In a very thin, 
suspended film with two free surfaces, the result is a 2D 
electroconvective instability, beyond which the fiuid flows 
in the plane of the film. 

The instability we consider is in principle possible in 
any isotropic, thin film of conducting fluid with a poten- 
tial difference applied across it. To realize this experi- 
mentally, we use thin suspended films of a liquid crystal 
in the smectic A phase. [l|, p^p^ 15 l^Jl^ These submi- 
cron films are newtonian and isotropic for fiows within 
the plane of the film, are immune to thickness variation, 
and have a very low electrical conductivity. The latter is 
required to keep the onset voltage within an experimen- 
tally reasonable range. The fact that convection occurrs 
in a smectic liquid crystaline material is irrelevant except 
insofar as this enforces its two-dimensionality. A simi- 
lar instability has been observed in suspended drops of 
other non-smectic fluids |Q , but these are much thicker 
and tend to experience large thickness variations when 
flowing, and so are much less 2D. It is also possible to 
convect other, more complex, smectic phases which are 
anisotropic in the film plane. We will not consider 

such anisotropy here. 



The system we consider consists of a thin film confined 
in the annular region between concentrically placed plane 
circular electrodes, as shown in Fig. |^. In the theory, 
both the film and the plate electrodes lie in the z = 
plane, and are treated as 2D surfaces. A voltage dif- 
ference is applied to the electrodes which define the in- 
ner and outer edges of the annulus. Experimentally, we 
find, in the absence of Couette shear, that the bifurcation 
from the quiescent and conducting state to the convecting 
state results in a flow pattern which consists of station- 
ary, counter-rotating pairs of vortices arranged around 
the circumference of the annulus. ||l[ Each pair of vor- 
tices has reflection symmetry and the overall pattern is 
spatially periodic and ID in the azimuthal direction. 

We alter the structure and symmetry of the quiescent 
base state by applying a Couette shear by rotating the 
inner electrode, as shown in Fig. |^. The shear flow leads 
to a net mean flow in the azimuthal direction, but unlike 
the previously studied cases of an imposed through-flow 
P, pl]p^ , the mean flow here is closed on itself. In 2D, the 
Couette flow by itself is stable |^ 21|, so that shearing 
alone does not introduce other instabilities. Experimen- 
tally, we find that the bifurcation from the Couette and 
conducting state to the convecting state leads to a ID 
pattern which travels azimuthally. |l| When this pattern 
is viewed in a co-rotating frame, that is, a frame in which 
the pattern is stationary, we observe that it is comprised 
of asymmetric counter-rotating vortex pairs. The vor- 
tices with circulation in the same sense as the rotation of 
the inner edge are narrower than those with the opposite 
sense. We will show below that these features of the non- 
linear flow patterns are already present in the unstable 
linear modes that we find from our stability analysis. 

An important effect of the Couette fiow is to suppress 
the onset of electroconvection. j|] The Couette base state 
is stabler against the radial electrical forcing than the un- 
sheared state. We will show that this additional stability 
is a feature of the shear and not of rotation. The central 
testable result of the linear stability analysis we present 
here is the degree of this suppression for various rates 
of shear. We directly compare this with our experimen- 
tal results obtained with annular smectic films and find 
satisfactory agreement. In addition, the linear stability 
analysis predicts the most unstable azimuthal mode and 
its rate of travel around the annulus. 

A crucial and somewhat subtle aspect of the phe- 
nomenon is the manner in which the electrical driving 
force couples with the fluid. Although the fiuid motion 
is confined to a plane, the form of the electric fields ex- 
terior to the plane, in the empty space above and below 
the film, are essential. These exterior fields are required 
for the self-consistent determination of the charge den- 
sities and thereby the electrical force. Since the electric 
potential at any point depends on the charge density ev- 
erywhere, the coupling of the two is inherently nonlocal. 
This is in sharp contrast to the case of RBC, in which 
the coupling of the density and buoyancy body force is 
entirely local. M 
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In Section ||, we develop our physical model of elec- 
troconvection in a conducting 2D fluid. This treatment 
is a generalization of that presented in Ref |^ . The base 
state solution of these electrohydrodynamic equations in 
an annular geometry is presented in Section [II. Fortu- 



nately, we find that the base state configuration of poten- 
tials and charges can be solved analytically in the annular 
geometry. The electrostatic configuration is independent 
of the Couette shear which forms the base flow. The 
shear is imposed by applying the appropriate boundary 
conditions on the flow velocity at the inner edge of the 
annulus. 

The linear stability equations of the sheared base state 
for non-axisymmetric perturbations are set up in Sec- 
tion 1^. Unlike for the base state itself, we find that 
no completely analytic solution exists of the electrostatic 
equations for the perturbed quantities. In the next two 
sections, we present two methods for solving the 3D, non- 
local coupling of the potentials and charges. In Section 
we make a local approximation to the electrostatic 
problem which is then used to close the perturbation 
equations. This yields approximate but useful analytic 
expressions for all of the field variables and the neutral 



stability boundary. In Section VI, we relax the local ap- 



proximation and use a numerical scheme to solve the 
electrostatic equations. These results are exact, up to 
the accuracy of the numerical scheme. We find that the 
local approximation does very well over most of the pa- 
rameter ranges of interest. For both solution schemes, 
we observe that shear suppresses the onset of convection 
and drives down the value of the azimuthal mode number 
of the convection pattern. 

In Section VII, we review our experimental apparatus 
for convecting and shearing liquid crystal films Q|, and 
compare the predictions of the previous sections with ex- 
periment. All of the material parameters of the fluid are 
known, or can be deduced from the data, so that the the- 
oretical results for the suppression of the onset of convec- 
tion can be compared to the data without adjus table pa- 
rameters. We find good agreement. In Section VII], we 



compare our results to what is known about other stan- 
dard instabilities , particularly RBC, under shear and/or 
rotation. Section [X is a brief summary and conclusion. 



II. THE GOVERNING EQUATIONS 

In this section, our physical model describing electro- 
convection in a suspended thin film is presented. Our 
analysis here is similar to that of Ref Unless oth- 
erwise stated, all operators, field variables, and material 
parameters are 2D quantities. 

The film is treated as a 2D conducting fluid in the 
z = plane, with areal density p, molecular viscosity 77, 
and conductivity a. We consider incompressible fluids, 
so that the 2D velocity field, u, is divergence free. 



Vu = 0. 



(2.1) 



The Navier-Stokes equation with an electrical body force. 



P 



^-f(u.V)u 



= -VP + 7?V^u + qE , 



(2.2) 



governs the fluid flow, where V, P, g, and E are the 
2D gradient operator, pressure field, surface charge den- 
sity, and electric field in the film plane, respectively. The 
term qE is the electric force acting on the surface charge 
density. The charge continuity equation 



dq 
'di 



= -V • {qvL + crE) , 



(2.3) 



takes into account the convective and conductive current 
densities, q\i and ctE respectively. Diffusion of charge in 
the plane of the film is neglected. 

Subscript three will be used to denote three- 
dimensional (3D) differential operators, material param- 
eters and field variables. The 3D electric potential \l/3 is 
governed by the 3D Laplace equation. 



0, 



(2.4) 



where V3 is the 3D gradient operator. The coupling of ^'3 
with the 2D charge density q is specified by requiring ^3 
to satisfy certain boundary conditions on the z = plane. 
The surface charge density q depends on the discontinuity 
in the z— derivative of 4*3 across the two surfaces of the 
film: 



q = -eo- 



= -2eo 



dz 



dz 



=0+ 



z=0+ 



z=0- 



(2.5) 



wher e ep is the permitivity of free space. If q is known, 
Eqs. |2.5| constitute Neumann conditions on ^'3 on the 
film, while Dirichlet conditions, described below, hold on 
the electrodes. If instead Dirichlet conditions are speci- 
fied on the film, Eqs. 2^ can be used to determine q. 

The 2D and 3D potentials are related via \1/ — ^'3|2=o. 
Equations 2.5 relate the surface charge density to the 
discontinuity in the z component of the 3D electric field 
E3 = — V35'3 across the film plane. On the other hand, 
the X and y components of E3 which form the 2D elec- 
tric field E = — V^*, are continuous across the film. This 
continuity is required by the usual matching conditions 
for electric fields across the surfaces of dielectrics. Note 
that it is the 2D quantity E and not E3 that appears in 
equations 2.2 and 2.3. One cannot simply use a Maxwell 



equation to eliminate the charge density in favor of the 
field because the 2D quantities in question are con- 
fined to a plane embedded in a 3D, otherwise empty, 
space and in gene ral V • E ^ g/eo. 



Equations 2.1- 2.5 



together with the appropriate 
boundary conditions on the electrodes, model our sys- 
tem. We consider the electrohydrodynamic limit where 
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magnetic fields and the resultant Lorentz forces are neg- 
ligible. One can also show that dielectric polarization 
effects are negligible in the limit of a thin film. |^ 

In our subsequent analysis, we define the stream func- 
tion by 



V X 



(2.6) 



where 

the pres sure field by applying the curl operator, Eqs. 



■z. Using Eq. 2.6, E — — V^* and eliminating 



and 2.3 



reduce to 



2.2 



• V 



(V X V X 



-r;V^(V X V X 0) + (Vg X V*) = , 



and 



dq ,^ 



Vq - fiV^* = 0. 



(2.7) 



(2.8) 



We render these equations dimensionless by rescaling 
the length, time, and electric potential by d, eod/cr, and 
y, respectively, where d and V are the cross-film width 
and potential difference. It follows that the stream func- 
tion and charge density are nondimensionali zed by crd/ert 
and e^V/d. Applying this rescahng to Eqs. 
and 2.5 gives 





Vdt 



(V X V X 0) + 7^^V^' X Vq^ 



^(^(Vx(?).vVvx Vx0), (2.9) 



dz 



(2.10) 

(2.11) 
(2.12) 



=0+ 



where, the dimensionless parameters 



7^ ; 



and P = 



£0% 



pad p^a^sd ' 



(2.13) 

are analogous to the Rayleigh and Prandtl numbers in 
the Rayleigh-Benard problem. Here, s is the thickness 
of the film and our 2D treatment assumes that s <^ d. 
The 2D material parameters are related to their three- 
dimensional counterparts by a = cr^s, rj = rj^s, and 
p — p^s. The control parameter TZ is proportional to 
the square of the applied voltage difference, but indepen- 
dent of the film width d. V is the ratio of the charge 
relaxation time scale in a film e^d/a^s to the viscous re- 
laxation time scalepad^/ 773. 

Equations 2.9-2.12, together with appropriate bound- 



ary conditions, describe electroconvection in a thin con- 
ducting film suspended in otherwise empty space, for any 
2D arrangement of the film and electrodes. 



III. THE BASE STATE 

We now specify the governing equations introduced in 
Section || to the case of an annular film, and solve them 
for the case of a general Couette shear fiow which forms 
the potentially unstable base state. 

We employ cylindrical coordinates (r, 0, z) as in Fig. 
The film is suspended between two circular electrodes 
which cover the remainder of the z = plane. The inner 
electrode has a radius r.i and is at potential 1 in our 
dimensionless units. The outer electrode, which occupies 
the 2 = plane for r > Tq, is at zero potential. The 
cross-film width is Tq — ri — 1 in dimensionless units and 
we define the radius ratio. 



so that 



1 



1 



(3.1) 



(3.2) 



Rotation of the inner electrode about the central r — 
axis produces a Couette shear in the base state. We 
will denote base state variables by the superscript zero. 
Under shear, the radial derivative of base state stream 
function is given by 



1 -a2 



1 



r(l - a)"" 



(3.3) 



where VL is the dimensionless angular rotation frequency 
of the inner electrode. If the fiuid is not sheared, the base 
state velocity field is zero and 0*^"^(r) = 0. The strength 
of the shear is described by a Reynolds number, with the 
velocity determined by the motion of the inner edge and 
the length by the film width. 



Tie = 



(3.4) 



where V is the Prandtl number given by Equation 2.13 



There is no loss of generality by treating only rotations 
of the inner electrode. Since our system is 2D, one can al- 
ways reduce independent rotations of both electrodes to 
this case by transforming to a rotating reference frame in 
which the outer edge is stationary. The transformation 



introduces a coriolis term into Eq. 2.2, which can simply 
be absorbed into the pressure. |^ As a result there are 
important differences in the stability of 2D and 3D sys- 
tems, and in particular for r igid rotation. We will return 
to this point in section VIII below. 



The base state potential '^'f^ (r, z) and charge density 
(7*^°) (r) are independent of the base state shear flow. They 
are determined by the electrostatic boundary value prob- 
lem given by Eqs. 2.11 and [2.12 , with Dirichlet boundary 
conditions on the z — plane 
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¥"\r) = *f^(r,0) = 



1 

_i_[ln(l - a) + ln{r)] 



< r < 
Ti < r < 
r >ro 

(3.5) 



The boundary condition for < r < To is found by 
treating the annular film as a 2D ohmic conductor subject 
to a dimensionless potential of 1 at the inner edge and 
at the outer, and requiring the continuity of the 2D 
current density. The logarithmic form follows from the 
cylindrical geometry. 



r(a = 0.8) 

4.25 4.5 4.75 



The Laplace equation 2.11 for the potential in the half- 
space z > is solved by the ansatz 



dk e'-'"'' Jo{kr)A{k) , 



(3.6) 



where Jq is the zeroth order Bessel function. Inversion 
of the above equation results in 



A{k) ^k dr r^f\r,Q)jQ{kr) . 
Jo 



(3.7) 



Hence, the base state surface charge density is given by 

D 

dC C*f (C,o)Jo(fcC) 

(3.8) 



Q(")(r) = 2/ dkk^Joikr) 
Jo 



where ^ is a dummy integration variable. Evaluating the 
integrals p§,pa , we find 



gW(r) 



2 
In a 



' 2'2' ' 



2' 2' 



(3.9) 



where F is the hypergeometric function 2F1. This func- 
tion is plotted in Fig. || for two values of a. As a ^ 1, 
9^°-' approaches the base state charge density for a later- 
ally unbounded rectangular film, which is odd-symmetric 
about the midline of the film. This symmetry helped sim- 
plify the analysis in that case. |^ However, for < a < 1, 
the annular base state charge density q^^^ is neither even 
nor odd about the midline, so the analysis here is more 
complicated. This deviation from odd symmetry is larger 
for smaller a. 

The surface charge density shown in Fig. |^ is "in- 
verted" in the sense that the positive charges lie close 
to the positive, inner electrode and negative charges are 
near the outer electrode. This unstable surface charge 
configuration gives rise to an electroconvective instabil- 
ity, much like the unstably stratified density configura- 
tion which produces Rayleigh-Benard convection. The 
divergences of q^*^^ at the edges of the film are a conse- 
quence of the idealized geometry in which the electrodes 
have zero thickness. 
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FIG. 2. Base state charge density g''''(r) at a = 0.5 and 
a — 0.8. g'"-* diverges at the edges of the film. 



This idealization leads to the boundary conditions in 
Eq. (3.5) , whose radial derivative is discontinuous at 
r = r, and r = r^,. This is reflected in Eq. 3.9 where 



F^i,i;l;l^ is indeterminate, so that g'^^ diverges at 
the edges of the film. 

IV. LINEAR STABILITY ANALYSIS 

In this section, we test the stability of the axisymmet- 
ric base state to non-axisymmetric perturbations. The 
perturbed quantities will be denoted by the superscript 
one, and are defined by 



c/,(r,0) =0(°)(r) (r,0), 
*(r,6l) = ^("^(r) + (r, 61), 



(4.1) 



*3(r,0,z) = *f (r,z) + *^^^(r,0,z) 



(1)^ 



Subs titution of the perturbed field variables into Eqs. |2.9- 
2.12 and retaining only the terms which are linear in the 
perturbed quantities yields 



V2- 



]_ d_ 
P\dt 



19(^(0) d 
r dr 89 



(V2</,(i)) 



~ ~ 00 00 



r \ dr 



(1) 



rV 36 Or \ dr^ r dr J 



(4.2) 
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where 



Q^2 



ld_ 

r dr 



1^ 









dr J 








(4.3) 




/A A\ 

(4.4) 


dz 


z=0+ 


(4.5) 


and 




" dz^ ■ 



The variables 4>'^^\ 5'^^^ and satisfy the following 
boundary conditions: 

0(i)(r„0) =5,0W(r„^?) = 

0(i'(^o,e) = 9.<^W(^o,e) -0, (4.6) 

*W(r„0) = *(i)(ro,0) =0, (4.7) 

*^^^(r,6l,z) ^ for z ^ ±00. (4.8) 



Equation 4.6 is a consequence of rigid boundary condi- 
tions on the fluid flow. The Dirichlet boundary condi- 
tions for the perturbed potential ^fg^-* on the z = plane 
arc 

{0 for < r < n 
*(i)(r,6l) for n<r<ro (4.9) 
for r > ro ■ 

We decompose the perturbations into products of the 
following form, 



/ 0(1) (r,0) \ 

*(i)(r, 6) 
^i'\r,9,z] 







y 






*™(r) 




V *3™(»-,2) / 



(4.10) 



where the azimuthal mode number m is an integer which 
corresponds to the number of vortex pairs in the pattern. 
The functions 0™, ^'m, and ^'3^ satisfy the same bound- 
ary conditions as the perturbations, Eqs. 4.6 - |4.9| . The 
growth rate 7 may be complex. 



Substitution of Eq. 4.10| into Eqs. |4.2[|4.5| gives 

77- 



7 - 



1 

V 
imTZ 



r 



ima 



rV 



TTi 

D,D- — I*, 



D.D--J- 



-D{D,D 

im 



(4.11) 



(4.12) 



D,D 



dz^ 



*3m = , 



= -2- 



3m 



dz 



(4.13) 



(4.14) 



z=0+ 



where D = dr and = D + 1/r. Henceforth, wc limit 
our discussion to a base state flow which is either qui- 
escent or Couette. In these cases D{D^D(t)'-°^) = and 
the right hand side of Eq. 4.11 is identically zero, It 
can be shown that in the limit a 
= 0, Eqs. 



4.11 



1, and for zero shear 
4.14 reduce to the linear stability 
equations for electroconvection in a laterally unbounded 
strip, dl 

We write the complex growth exponent 7 as 7''+77'. In 
order to find the conditions for marginal stability, we set 
the real part 7'' = 0. Our task is then to solve Eqs. [4.11[ - 
4.14 for a given a, V and 72.e, by determining consistent 
values of TZ and 7* for each 777,. The rate of azimuthal 
travel of each marginally stable mode around the annu- 
lus is 7*/77i. The marginal stability boundary, which is 
defined only at discrete 771, has a minimum at the critical 
values 777c, T^ci while the critical mode travels at ^^jnic- 

We employ the following basic expansions: 



^m{r) = ^ ^„ (7-) , 


(4.15) 


*m(r) = ^ A„*,„;„(r) , 

n 


(4.16) 


*3m(f', -2) = X! ^«*3m;ri(»', z) , 

n 


(4.17) 


Qmif^) — ^ ^ ^nQm;n(^) : 


(4.18) 



where the An are amplitudes. The expansion eigenfunc- 
tions 0m;„(r), 5'„i;„(r), and \E'3m;„(r, z) satisfy the fol- 
lowing boundary conditions: 

<i>m;n{ro) = dr4>r,i-n{r o) = , (4-19) 

*m;„(r-,;) = «'„.;„(ro) = 0, (4.20) 

^'3m;n(?', z) ^ for Z ^ ±00 , (4.21) 

( for < r < r, 

*3m;n(»-, Z = 0) = <^ *,„;„(r) for < T < To . (4.22) 

I for r > To 

The functions 4>m;nir), which satisfy the 'rigid' boundary 
conditions given in Eq. can be identified with the 

Chandrasekhar cylinder functions 0, 

~^Cfyin-^m{Pm7i^) ~^ ^jnn-^jniPjnn^) ■ (4.23) 



The boundary conditions Eq. 4.2C, imply that the 2D po- 
tential expansion function '^m-n can be further expanded 
in a series of functions of the form. 
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(4.24) 



and tprn-pi along with their associ- 
" ' ' are 



The functions Cm-.n 

ated constants Prnn ; ^mn ; Crnn 7 ^mn 7 \mp 7 and bj^p 

described in detail in the Appendix. 

In the next two sections, we use these expansions to 
solve Eqs. |4.1l| - [4.14 , first approximately and then in gen- 
eral. 



V. LINEAR STABILITY: LOCAL 
APPROXIMATION 



The main barrier to solving Eqs. 



4.11^4.141 lies in 
the difficu lt nonlocal coupling between and qm in 
Eqs. 4.11 and 4.12 whi ch is requ ired by the 3D elec- 
trostatic equations 4.13| and 4.14 . In this section, we 
circumvent this problem by making an approximation. 
This yields expressions which give some insight into the 
general linear stability problem. 

To make the approximation, we replace the 3D 



relation: 



Eqs. 4.13 and 4.14 with the following simple 2D closure 



(5.1) 



In the above expression, is a closure factor which is 
independent of r, and is to be specified. As we show 
below, a consequence of this approximation is that the 
charge density and the 2D potential are related pointwise, 
or locally. 

This approximation is motivated by following physical 
reasoning. If instead of an annular film, one considers 
an annular column, with a height much larger than its 
width, then there is a straightforward Poisson relation 
between a bulk charge density and the 3D potential in- 
side the column. If the 3D potential is independent of 
z, and is equal to 5'm(r)e™^, there are no free surfaces 



to consider and one has in place of Eqs. 4.13 and 4.14 
the relation (^D^D — m^/r^^'I'„i = —qm- If one now hy- 
pothesizes that the charge density retains its radial profile 
when the bulk is 'squeezed down' to a 2D ffim, then one 
must include only an r independent scaling factor, /„, 
as in Eq. 0|. 



We ma ke us e of 4.16 and further expand each ^E*^ 
using Eq. |4.24| so that 



(5.2) 



where the Bm-np are constants. It then follows from 



Eqs. M and Eq. ^ that 



qm — ^ ^ ^ ^ Bm;np fm;p)Cmp ^m\p — 

^ ^ ^ ^ Bm;npqm;p 7 (^■*^) 



n p 




FIG. 3. The charge density perturbation qa-iir) at 
a = 0.50, as determined by the local approximation of Sec- 
tion ^ (dashed line) and the nonlocal solution of Section ^ 
(solid line). The nonlocal result diverges to —00 at the edges 
of the film. 



where 



qm]p — fm-p/Cmp '^m-^p 7 



(5.4) 



which demonstrates the pointwise, local relation between 
the potential on the film and its surface charge density 
in this approximation. 

A choice for /m;p can be made by considering the 
empty upper half space z > with homogeneous bound- 
ary conditions at infinity and Dirichlct boundary con- 
ditions such that the potential is equal to '4'm;pir)e'^™'^ 
everywhere on the z — plane. This boundary condi- 
tion is sm ooth, unlike the piecewise smooth conditions in 
Eq. 4.22, which respect the sharp edges of the annulus. 



Then the 3D potential for z > satisfies the equation. 



D^D- 



(5.5) 



Using the eigenvalue relation Eq. A9, it follows that 



is thus 



Xmp- The corresponding surface charge density 



oz 



'^Xmpi'm-p- (5.6) 



z=0+ 



Comparing Eq. 5.6 with Eq. 5.4 leads to the following 
choice for the closure factor: 



m;p 



Xmp 



(5.7) 



n p 



For the remainder of this section, we adopt this method of 
closing the equations. Fig. || shows a plot of the approx- 
imate charge density, qe-i at a = 0.5, corresponding to a 
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potential i/jq-i on the film. It is compared to a more ac- 
curate numerical solution discussed in Section VI, below. 



As might be expected, the approximation is accurate ex- 
cept close to the edges of the film. 

We can now solve the remaining 2D e quat ions using the 
expansions in Eqs. 4.15 , 4.23 , 5.2 and |5.3| . We consider 
the simplest case: a single expansion mode, so that Ai — 
1 and An — for n > 1. Similarly, we truncate the 
expansions of the potential and charge density at p — 1 
E for p > 1. The expansion coefficient 



so 


that Bm 




= Bjn-l 


5.2 


and [5^ 



— Bjnai is complex. Thus, using Eq. 5.7, Eqs. 4.15 



4^m — ; 1 1 



(5.8) 



Substitution of Eq. 5.8 into Eq. 4.12 , with 7 = 17' results 
in an equation, the real and imaginary parts of which are 

B^nXmi^n.:! = 2 f f - a^E^^) Bl^i,^,^ , (5.9) 



BlnXml V'm;l+2(7' j B^^jXmlV'm;! 

+ (^)C™;1 = 0, (5.10) 

respectively. Eliminating B^^ from the above pair of 
equations results in 



4 / , mD 
1 + 2 



m;l 

C„,;l = 0. (5.11) 



Multiplying by ipm-,!, and integrating to form inner prod- 
ucts. 



(...) 



rdr . 



we solve for the expansion constant. 



(5.12) 



(5.13) 



where -hftp^.-i is a normalization factor given in the Ap- 
pendix and the matrix element 



(5.14) 



A similar projection of Eq. 5.9 and some simplification 
results in 



Bl 



2^r« 
Xml 



D 



(5.15) 



Equations 5.13| - ^.15 determine the expansion coefficients 
of the potential and charge density for a giv en stream 
fimction. Substitution of Eq. 5.8 into Eq. 4.11, with 7 = 
^7* and a Couette shear gives 



ft ^ri _ ■ Pml 



I (Vm;i — Um-i) 



rriTZ 



— (b;„ + zs;j(2 



-i5g(°))V'™;i j = 0, (5.16) 

where lAm-^i and Vm;i are defined in the Appendix. Pro- 
jecting the above equation with C,„;i allows us to solve 
for U: 



M^Bln - Bin) 2 XmlJm - K 



(5.17) 



where 



F — 


(c...."f°'(^ 






Jm — 









(5.18) 
(5.19) 

(5.20) 

(5.21) 



The normalization factor Mc^ -l is given in the Appendix. 

To determine the linear stability boundary for a given 
radius ratio a, Pr andt l number V and Reynolds number 
T^e, we solve Eq. 5.17 for a sequence of azimuthal mode 



numbers m, using Mathematica for all integrations. We 
first consider some special cases. 
For zero shear, 



72. (a, m) 



= 7* = and Eq. |5.17| reduces to 

/3ml''Xml^-^C™;iA/'^™;pi 



m2K„i 2 Xml Jm - K„ 



(5.22) 



Note that Eq. 5.22 is independent of the Prandtl number 
and is always real. Fig. ^ shows the resulting stability 
boundary for a = 0.5. Except for being discretized in in- 
teger values of m by the annular geometry, the boundary 
is not much different than that of the infinite rectangular 
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200 - 




m 

FIG. 4. The marginal stability boundary TZ vs. m for zero 
applied shear and a = 0.5. Open symbols show the result 
from the local approximation (Section ^ , solid symbols show 
the full nonlocal stability calculation (Section [Vl|) , with three 
expansion modes. The dashed line indicates the continuous 
marginal stability curve for the unbounded rectangular ge- 
ometry with plate electrodes (see Ref. ^). The azimuthal 
mode number m is equivalent to a dimensionless wavenumber 
K = ra/f (upper scale), where f is the mean radius (ri + ro)/2. 

case which is defined at a continuum of wavevectors 
K. The discrete curve approaches the continuum one in 
the hmit a ^ 1. Also shown in Fig. |^ is the more accu- 
rate numerical solution presented in Section V], below. 
The minima of this curve define the critical parameters 
TO°(a) = mc{a,ne = 0) and 7^2(a) = 7^c(a,7^e = 0). 
Some values for various a a re co llected in Table |[ 

For non-zero shear, Eq. ^.17] simplifies somewhat in 
the limit that the Prandtl number V oo. A word of 
caution is perhaps appropriate here. As "P — > oo, we 
require ^ oo such that fi/P and therefore TZe (see 
Eq. are fi nite . In this limit, with the proviso that 
TZ be real, Eq. 5.17 becomes 



72. (a, TZe, m, 7') 



/3ml'^Xml^A/c„,.iA/'^„,.yi 



(5.23) 



The only shear dependence in Eq. 5.23 occurs thr ough 
the matrix element L^, which is defined by Eq. ^.14[ 
For zero shear, L™ = Km, while for non-zero shear, 
Lm is bounded above by K„j. Hence, TZ{a,TZe,m,^^) 
is bounded below by TZ{a,TZe ~ 0,m,7* — 0), the zero 
shear value. Thus, one expects suppression of the onset of 
convection for non-zero shear for every non-axisymmetric 
mode TO. This feature is also present for finite V. 



radius 


local 


nonlocal 


nonlocal 




ratio 


(n = 




(n = 


-1,1 = 20) 


(n = 


= 3,Z = 


20) 


a 




n'l 






m'i 






0.33 


4 


102.58 


4 


91.62 


4 


4.03 


82.15 


0.56 


7 


91.25 


7 


92.47 


7 


3.95 


89.64 


0.6446 


9 


90.04 


10 


93.00 


10 


4.33 


91.12 


0.80 


18 


88.84 


19 


93.59 


19 


4.22 


93.10 



TABLE I. Critical parameters for zero shear as determined 
by the local approximation (Section ^ and the full nonlocal 
(Section VI) linear stability analysis. The integer n (p or Z) is 
the number of modes (expansion functions) used in the series 
representation of the field variables. The critical wavevector 
for electroconvection in a laterally unbounded geometry is 
Kc = 4.223 (see Ref. Q). The ratio m!i/f — > Kc as a — > 1, 



where r is the mean radius (r^ + ro)/2. 



Prandtl 


Rotation 


local 




nonlocal 






rate 


(n = 




1) 


{n-- 


^1,1 = 20) 


V 


Q. 




7^c 


7c 


nic 


7^c 


Yc 


100 


5.5135 


8 


319.50 


-14.949 


8 


282.61 


-16.383 


10 


0.5513 


9 


99.65 


-1.655 


9 


98.52 


-1.861 


1 


0.0551 


9 


90.16 


-0.166 


10 


93.06 


-0.205 


0.1 


0.0055 


9 


90.05 


-0.017 


10 


93.00 


-0.020 


0.01 


0.0006 


9 


90.04 


-0.002 


10 


93.00 


-0.002 



TABLE II. Critical Parameters for a range of Prandtl num- 
bers V at radius ratio a = 0.6446 and Reynolds number 
TZe = 0.1. The integer n (p or /) is the number of modes 
(expansion functions) used in the series representation of the 
field variables. 



Prandtl 



P 



100 
10 
1 

0.1 
0.01 



Reynolds 



Re 



0.0181 
0.1814 
1.8137 
18.1373 
181.3731 



local 

(n = l,p = l) 



115.20 
115.29 
115.27 
115.24 
115.16 



-2.995 
-3.002 
-3.007 
-3.039 
-3.086 



nonlocal 

{n = l,l = 20) 



7^c 



109.03 
109.02 
108.99 
109.14 
109.31 



-3.372 
-3.365 
-3.307 
-3.165 
-3.108 



TABLE III. Critical Parameters for a range of Prandtl 
numbers V at radius ratio a = 0.6446 and rotation rate 
57 = 1.0. The integer n (p or I) is the number of modes 
(expansion functions) used in the series representation of the 
field variables. 
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For arbitrary a, V and TZe, one can solve Eq. 5.17 for 
various to by a ID search procedure. At each to, 7' is 
varied to find a TZ that is real and a minimum. Example 
neutral curves are shown in Fig. ^ for a = 0.8, P = 10 
and several TZe. The suppression of convective onset is 
evident, as well as a tendency for the critical mode num- 
ber TOc to decrease with TZe. We find that TZ{TZe) is 
a monotonically increasing function of TZe. The depen- 
dence of 7^ on T' is treated in two ways. We will first fix 
TZe, change V and consequently let vary. In the sec- 
ond protocol we will fix f2, change V and consequently 
let TZe vary. For fixed TZe, the P dependence of these 
results is not strong, except at large T'. Because TZe is 
proportional to V^^ in Eq. 3.4, this limit corresponds to 
large values of fl and a large suppression effect. Some re- 
sults for various T' are summarized in Table ^ For fixed 
51, the P dependence is very weak at all V investigated. 



Results are summarized in Table III 



VI. LINEAR STABILITY: EXACT NONLOCAL 
SOLUTION 

We now relax the local approximation described in the 
previous section, and properly treat the full problem. 
The electrostatic Eqs. 03-1TJ, in which the charge den- 
sity and electrostatic potential are related nonlocally, are 
solved numerically. As well as being more physically cor- 
rect, this allows us to determine the accuracy of the local 
approximation. 

The first step in the present method is to find the ap- 
propriate expansion functions of the field variables. Sub- 
stitution of Eqs. |4.15| - [1.18| into Eqs. ^1.12[j4.14| yields equa- 
tions which can be solved for A„, (j)m,-n, ^'m;«, ^Zm-n, 
and Qm; n- The stream function (r) — Cm;nir) as in 
Eq. 4.23 , The potentials ^'m;n and ^'3„i;n and the charge 
density Qm-n are further expanded as 



^m-nir) = ^ Vm-nlfpra-lir) , 
I 

I 

Qm-n{r) = ^ nm-niqrn;l{r) , 



(6.1) 
(6.2) 
(6.3) 



where ipm-j is given by Eq. 4.24 and Vm;ni are complex 
coefficients. The functions tpsm-i and qm;i are computed 
as follows. After substituting Eqs. 4.17 and 6.2 into 



Eq. 4.13, the resulting equation 



777,^ 0'^ 



(6.4) 



is solved numerically on a finite 2D grid by an over- 
relaxation algorithm p3| for the functions il' sm-.ii ' r, z) 



4.21-4.22 



subject to the boundary conditions as in Eqs 
with ip3m-iir,0) = ipm :i(r) for r,, < r < r„ and othcr- 
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FIG. 5. Marginal stability boundaries under shear, showing 
TZ a.s & function of m for a = 0.80 and P = 10, using the local 
approximation (Section Shown are the boundaries for 

Reynolds numbers TZe =0.0 (circles), 0.25 (filled circles) and 
0.75 (squares). Arrows indicate the critical mode. 



(6.5) 



where the differentiation is performed numerically. Fig. ^ 
shows a plot of the charge density gga at a = 0.5, and 
compares the numerical result to the approximate one 
given by Eqs. 5.4 and 5.7 from the previous section. 



The approximate solution does not contain the diver- 
gences which occur near the film's edges, due to the sharp 
changes in the derivative of the potential. These are a 
feature of our 2D model, which treats the film and elec- 
trodes as having zero thickness. 



We next substitute Equations S.l and 6.3 into Eq. 4.12 
and use 



Cn 



I r 




■ I )Vni;l ■ 



(6.6) 



wise. Then Eqs. 4.14, 4.17, 4.18, and 



4 give 



The resulting equation is projected against 'ipm;k to ob- 
tain a matrix expression that can be solved numerically 
for the complex coefficients Vm-ni- 

Finall y, sub stituting the various expansions in ipm-i 
into Eq. 4.11 and taking the inner product with Cm-p, 
yields a set of linear homogeneous equations for the con- 
stants An- We write this set as the matrix equation 
J2n ^nTp„ = 0. For a nontrivial solution, the compati- 
bility condition is 
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Real{det [T]) = Imag(dei [T]) = 0, 
with the elements of the matrix T given by 

2 



(6.7) 



pn 



m 



r 



The first two inner products of Eq. 3.8 can be simplified 
using Eqs. A_l and |A8| and the matrix elements F^i;^^ are 



■ m;pn 



E 



a 



m\p 



iD¥'>^)q^,i-{Dr>)i,„,a 



(0)) 



(6.9) 



The real values of TZ and 7* which satisfy Eqs. 6.7 and B 



at each m define the neutral stability boundary TZ = 
TZ{a,V,7le,m,Y)- The critical parameters mc, TZc and 
7* are obtained when TZ is minimized. 

The numerical over-relaxation calculation used to solve 



Eq. 6.4 involved a grid spacing such that there was a min- 



imum of 160 points across the width of the film. For the 
purposes of integratio n, th e discrete values of Qm-i found 

wer e Ch ebyshev interpolated. 
-|6.3| and |6.6| were calculated up to 
p — were employed in the 



numerically from Eq. 6.5 



The series in Eqs. 6.1 
^ = 20. Three modes (n 
compatibility conditions Eq. 6.7 when the shear was zero. 
This was reduced to one mode {n = p = 1) when the 
shear was non-zero. All r-integrations were performed by 
the Romberg method. ||§ We solved Real{det [T]) = 
(Eq. |6.7D for TZ for a given trial 7* and associated coeffi- 
cients Vjn-ni of the field variables. The TZ and trial Vm-ni 
were employed in the search for the 7* which satisfied 
lmag{det [T]) = 0. The new value of 7* was then used 
to find new coefficients Vm-ni- The iterative cycle was 
continued until the parameters and field variables were 
converged. 

The marginal stability boundary for a — 0.5 and zero 
shear is shown in Fig. using both the local approxi- 
mations and the full nonlocal solution. Neither solution 
is qualitatively much different from the infinite rectan- 
gular case. 1^ For the nonlocal solution, we find that as 
a ^ 1 , mP^ increases such that m° / f approaches the cor- 
rect limiting value, Kc = 4.223. Here, f = (r^ +ro)/2 
is the midline radius of the annulus and Kc is the crit- 
ical wavenumber for an infinite rectangular film in the 
"plate" electrode geometry. ||^ Table || shows a sample 
set of critical parameters as determined by both solution 
schemes. There is generally good agreement between the 
two methods, except at small a, where more expansion 
modes are needed. 

Fig. ^ shows the zero shear critical control parameter 
TZ^ as a function of the radius ratio a. There are discrete 
values of a where 72.^ (a, TOc) = TZ'^ (a, rric + 1), so that 




FIG. 6. The critical control parameter for zero shear TZ^ 
versus radius ratio a using the nonlocal stability analysis with 
three expansion modes. Special radius ratios, at which two 
adjacent modes m and m + 1 are simultaneously marginally 
unstable, occur at the cusps of the curve. Between cusps, the 
critical mode number remains constant. 

two adjacent modes become unstable simultaneously at 
a co-dimension two point. These points occur at the 
cusps in Fig. ^, while between the cusps a single value of 
m is critical. It is interesting to note that the trend in 
TZ'^ is increasing overall, as function of a. This is oppo- 
site to what is found for radially driven Rayleigh-Benard 
convection in a rotating annulus. |^ We attribute this 
difference as a manifestation of the differences between 
the nature of the charge density 'inversion' and thereby 
electrical forcing with the bouyancy inversion of RBC. 
As a ^ 1, the co-dimension two points become closely 
spaced and the value of TZ^ approaches a limiting value. 



For the case of non-zero shear, the nonlocal analysis 
produces neutral curves which resemble those shown in 
Fig. ||. For fixed a and 'P, the critical control parameter 
TZc increases with TZe, which indicates a suppression of 
convection as is observed in experiments. This pre- 
diction is compared to experimental data in Section VII 
below. 



The P dependence of TZc for TZe 7^ is rather 
weak at small 7^. This is outlined in Tables || and III 



Note that for fixed TZe, large P implies large and thus 
greater stabilization. For fixed Q, however, the stability 
propoerties are insensitive to T' and therefore to TZe. 

The critical mode number rric decreases with increasing 
shear. This effect is shown in Fig. |^, which displays 
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FIG. 7. Ratio of imaginary part of the growth rate to crit- 
ical mode number, Yc/^c versus Reynolds number TZe for 
a =0.56, 0.6446, and 0.80 and V = 123 using the nonlocal 
linear stability calculation with one expansion mode. The 
breaks in the curves for each a show the intervals over which 
the critical azimuthal mode number iric has the value indi- 
cated. 

the angular traveling rate of the critical mode, 7*/mc, 
as a function of TZe for several a. For > 0, 7* < 0, 
which indicates that the critical mode travels around the 
annulus in the same sense as the inner electrode, jl/rric 
is a very nearly linearly decreasing function of TZe, with 
very small discontinuities at points where the critical az- 
imuthal mode number rric changes, as shown in Fig. |^. 
Each of these discontinuities, which are too small to re- 
solve on the scale of Fig. |^, is a co-dimension two point, 
where two adjacent m modes with very slightly differ- 
ent traveling rates are simultaneously unstable at onset. 
These special points are also slightly P dependent, as 
well as being a dependent in a manner similar to the 
zero shear case discussed above. 

Although linear analysis cannot provide the magnitude 
of the fields above onset, it is nevertheless interesting 
to examine the spatial structure of the linearly unsta- 
ble modes. Figure || displays the velocity vector field 
of a critical mode for a = 0.56, plotted with an arbi- 
trary amplitude. In Fig. ^(a), we show the stationary 
vortex pattern at TZe — 0. Here = 7, so 7 symmetric 
vortex pairs are arranged around the annulus. This pat- 
tern is purely non-axisymmetric or "columnar" . It is an 
exact solution to the governing equations and has been 
observed experimentally. As is discussed in Section VIII 



below, we would expect exactly the same pattern under 
the same conditions for any rate of rigid rotation. In 
Fig. ||(b) , we show the typical flow pattern for 



(b) 



(c) 




FIG. 8. The velocity vector field of arbitrary amplitude 
for a = 0.56 with (a) TZe = 0, (b) TZe > as viewed in 
the laboratory frame, and (c) TZe > as seen in the frame 
co-rotating at -yl/mc in which the pattern is stationary. 

a large TZe > 0, as viewed in the laboratory frame. The 
periodicity of the pattern is reduced {i.e. iric < m'^), and 
the traveling pattern appears as a meandering wave in 
the laboratory frame. Fig. ||(c) shows the same TZe as in 
(b), but as seen in the frame in which the pattern is sta- 
tionary. This frame rotates in the same sense as the inner 
electrode but with an angular speed of jl/rric, which is 
less than fi. In this frame, each vortex pair consists of a 
larger and a smaller member. The same qualitative fea- 
tures are seen in sheared convection patterns just above 
onset in the present experiments and those of Ref . M . 



VII. EXPERIMENT 

In this section, we discuss our experimental investiga- 
tion of the linear stability of 2D fluids to electroconvec- 
tion and Couette shear. The apparatus we describe is 
only slightly modified from that previously reported in 
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Ref Our objective here is to experimentally test one 
of the direct predictions of the linear theory, namely the 
suppression of onset by shear. 

A smectic A liquid crystal (octylcyanobiphenyl, or 
8CB) was employed to realize a 2D fluid. In the smectic A 
phase, long organic molecules arrange themselves into 
layers with their long axes normal to the layer. Within 
the layer planes, smectic A behaves as a 2D isotropic 
fluid, and can be drawn into robust submicron films 
which are integer numbers of layers thick and uniform 
over their whole area. [Other smectic phases form 

2D anisotropic fluid films^ 3l| , p5| , p6t Unlike soap films, 
smectics can flow in the film plane without thickness 
change and are not susceptible to evaporation. Smectics 
generally have low electrical conductivity due to residual 
ionic impurities. | |3l| , |3^ This combination of properties 
makes smectic films uniquely suited for studying elec- 
troconvection. While our experiments used liquid crys- 
tal films, the electroconvective instability itself is generic 
and not specific to anisotropic fluids; a similar instability 
is found in thin drops or puddles of ordinary isotropic 
fluids. [ p4| , p3|j3^ ] The anisotropy introduced by the smec- 
tic layering simply serves to constrain the film to have a 
constant, very small, thickness. It behaves isotropically 
for flows in the plane of interest. 

Our experimental apparatus consisted of two concen- 
trically placed stainless steel electrodes which supported 
a suspended film. The inner electrode was a disk of radius 
3.61 ±0.01 mm and the outer electrode was a plate with a 
circular hole of radius 6.40 ± 0.01 mm, giving an annular 
film with radius ratio a = 0.564 ± 0.002. Each molecu- 
lar layer of smectic A 8CB is 3.16 nm thick. The inner 
electrode could be rotated about its axis by a high preci- 
sion stepper motor at angular frequencies up to 4 rad/s, 
while the outer was fixed. All experiments were per- 
formed at 23 ± 1°C. The air surrounding the film was 
pumped down to an ambient pressure of 0.5-5 torr. The 
low pressure ensures that air drag on the film will be in- 
significant. The system was electrically shielded by 
a Faraday cage and the outer electrode and shield were 
maintained at ground potential while the voltage V of the 
inner electrode was varied above ground. The resulting 
current through the film was measured using a picoam- 
meter. Except for a few qualitative runs to observe the 
vortex patterns, we did not visualize the flow, as this re- 
quired suspending particles in the film, which was not 
necessary for our present purposes. The film behaved as 
an ohmic conductor in the conduction regime below the 
onset of convection. The critical voltage Vc, at which the 
onset of convection occurred, was located from the sharp 
increase in current. In all runs, the shear rate was held 
constant while the voltage was slowly increased and de- 
creased through onset. For some rates of shear, the onset 
of convection was hysteretic. [|l| In these cases, we took 
Vc to be the transition from conduction to convection as 
V is increased. 

In order to interpret the experimentally measured 
quantities, we need to accurately express them in nondi- 



mensional form. According to Eq. 2.13, the critical con- 



trol parameter TZc is related to the critical voltage Vc by 



Uc 



f'^V^ 



(7.1) 



The size of the experimental uncertainties on a^, s, 



especially, rya prevented us from directly using Eq. 7.1 



and 
to 



scale the voltages. Rather, we used the zero shear mea- 
surements of Vc to scale the data for non-zero shear. A 
relative measure of the suppression of the onset of con- 
vection is 



e (7^e) = 



7^c(7^e) 



7^o 



- 1 = 



VciUe) 



1 2 



- 1 



(7.2) 



The first part of Eq. 7.2 was applied to the theoretical 



results, while the second was applied to the experimental 
data. The latter scaling does not require a knowledge of 
the material parameters. In order to obtain the remain- 
ing dimensionless quantities, TZe and V, estimates of (73, 
s, and 773 are required. However, we will show below that 
our final results are independent of s and (73 and quite 
insensitive to the value of 773. 

It is straightforward to measure a quantity propor- 
tional to the product (73 s from the data. For V < Vc, 
the film behaves ohmically and hence the slope of the 
current-voltage data is the film conductance c, given by 



27rCT3S 
ln(l/a) 



(7.3) 



For a given current-voltage sweep, c can be measured 
to within a few percent. This precision is limited by 
the tendency of the conductance of the material to drift 
over the course of the sweep. |3^ Using many measured 
values of c, we found an average value of , which was 
consistent with all the data. The film thickness s can 
be independently deduced by examining the color of the 
film under reflected white light. p6t By measuring c and 
s for a number of films with various a, we found the drift- 
averaged conductivity to be 173 = 1.3±0.2 X 10^^(rim)~^. 
We used this value to infer a consistent value of 773, via 
the procedure desc ribed below. 

Eqs. 7T and 7^ can be combined to give 



27reo 



ln(l/a)i 



(7.4) 



To find 773 , we plotted the left hand side of Eq. [7.4| against 
measurements of c, using measured values of V^ and the 
theoretical value of TZc- We used data from a large num- 
ber of films with s < 60 smectic layers for two values of 
a, 0.56 and 0.64. A on e para meter fit to find the slope of 
this plot established a/ 773 /173 . Using the average value of 
(73 given above, we found 773 = 0.16 ± 0.06 kg/ms. This 
method has the advantage that it does not require an 
independent measurement of s for every film. 
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FIG. 9. Comparison of the experimental and theoreti- 
cal relative suppression e versus Reynolds number TZe for 
a — 0.56. The dashed (solid) line is the result of the local 
(nonlocal) method. 

We then employed the above value of 773 to calculate 
TZe which is expressed in terms of dimensional parameters 
by 



TZe = p3LLiri{ro - ri)/ri3 



(7.5) 



where uj is the measured angular frequency of the inner 
electrode in rad/s. The 3D density of the liquid crystal 
1.0 X 10^ kg/m^. Similarly, the Prandtl number T^ 



1 37 



is given in terms of measured quantities by 



P = 



27reo?73 



1 



Psiro - ri)ln(l/Q;) 



(7.6) 



Fro m t he definition, Eq. 2.13, T' depends on s. In 



Eq. 7.6, this dependence is contained in c, which can be 
accurately measured for each individual current-voltage 
sweep. Thus, for films of various constant thicknesses, P 
experiences a small drift from sweep to sweep, due to the 
drift in c. For the data we report here, P was always in 
the range 52 < P < 75. If c is treated as an indepen- 
dent, experimentally determined variable, then Eq. 7.6 is 
independent of s and (T3. 

We can now compare our experimental results and the 
predictions of the linear theory for non-zero shear. Fig. ^ 
shows the theoretical curves from the local and nonlocal 
calculations and the experimental data for the relative 
suppression e, as a function of TZe. The theory curves 
are plotted for a fixed P = 63.5. There are no adjustable 
parameters. We find good agreement between the exper- 
imental data and the linear theory for TZe < 0.22. The lo- 
cal approximation slightly overestimates the suppression, 
while the exact nonlocal result slightly underestimates 



it. Recall that the local approximation neglects the di- 
vergence in the perturbed charge density at the edges of 
the film, but uses the full base state solution, which has 
similar divergences in our 2D model. The nonlocal calcu- 
lation keeps the divergences at both orders. The fact that 
these two treatments bracket the data suggests that the 
divergences are rounded off by the finite thickness of the 
film and electrodes in the actual experimental geometry. 

The good agreement shown in Fig. ^ is essentially in- 
dependent of the value of 773. Since the 773 dependence in 
the TZe scaling of both the theor y (v ia V in Eq. 3.4) and 



the experiment (according to Eq. 7.5) are proportional to 
1/773, any change in 773 multiples both by the same factor. 
This simply results in a rescaling of the TZe axis in Fig. ||, 
with no change in the quality of the agreement. 



VIII. DISCUSSION 

It is interesting to compare our results for 2D electro- 
convection under shear with those for other, more famil- 
iar instabilities under shear or rotation. 

In a theoretical study, Agrait and Castellanos |l^ con- 
sidered electroconvection due to a radial field in a 3D 
concentric cylinder geometry with charge injection on ei- 
ther cylinder. Both cylinders were permitted to rotate to 
produce a general Couette shear. Their result was that 
shearing enhanced the instability, leading to a 3D flow 
that resembles TVF. This can easily be understood as a 
consequence of the fact that in 3D, the base state Couette 
shear itself is unstable. In our 2D system, the shear flow 
is not only stable by itself |^,^, but, as we have shown, 
has a stabilizing effect on the electroconvection. This is 
one of several examples of how the 2D nature of our flow 
leads to an unusual stability, which is often difficult or 
impossible to realize experimentally in other systems. 

The added stability in our system is a consequence 
of the shear and not of rotation. Under rigid rotation, 
where the inner and outer electrodes are co-rotating, one 
can transform to rotating co-ordinates in which the elec- 
trodes are stationary. This transformation introduces a 
coriolis term —2Uz x v = —2ilWtlj in Eqn. 2.2 which 
may be absorbed into the pressure gradient term VP 
and eliminated. |^ Thus, in a purely 2D system, rigid 
rotation and the non-rotating, unsheared case have iden- 
tical stability. It also follows, since the transformation is 
general and the unsheared bifurcation is stationary, that 
the resulting nonlinear vortex pattern above onset must 
be stationary in the co-rotating frame. 

This lack of dependence under rigid rotation may be 
contrasted with a large class of 3D and quasi-2D rotat- 
ing Rayleigh-Benard systems where rotation 
produces added stability but the absence of strictly 2D 
flow results in time-dependence (precession) of the con- 
vection pattern in the co-rotating frame. Chandrasekhar 
treated the classic problem of the linear stability of 
RBC in a laterally unbounded layer rotating about its 
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normal. The case of a laterally bounded cylindrical layer 
has received much recent interest theoretically and 
has also been the subject of a precise experimental study 

A concise summary of some of the work on rotat- 
ing RBC can be found in Ref. Due to the simi- 
lar symmetries in these scenarios, some results are com- 
mon to most rotating RBC systems. Since, these 
systems are three-dimensional (3D), the solutions they 
support can generally be classified into axisymmetric, 
non-axisymmetric ("columnar") or mixed (combinations 
of axisymmetric and non-axisymmetric) solutions. The 
principal result is that for 3D mixed solutions, the onset 
bifurcation is no longer steady as it is for non-rotating 
RBC. This leads to a flow pattern that precesses in the 
co-rotating frame. The onset of these solutions is however 
suppressed, raising the critical Rayleigh number above 
its non-rotating value. [^j9|,[Tot Purely non-axisymmetric 
or columnar solutions are strictly two-dimensional (2D). 
When they do occur, they do not precess but their onset 
occurs at the same critical Rayleigh number as in the ab- 
sence of rotation. There have also been theoretical 
studies of the interesting but nonetheless experimentally 
unrealizable situation of 2D RBC in a rotating annu- 
lar geometry with purely radial gravity and heating 
These studies found similar columnar solutions. In fact, 
purely columnar solutions ("Taylor columns" |^) have 
yet to be observed in any rotating RBC experiment since 
the boundary conditions at the top and botton of the 
cylinder must be stress free |6|,|l0|, a requirement that 
cannot be attained in terrestrial RBC experiments. In 
contrast, two-dimensionality, stress free end boundary 
conditions and radial driving forces all arise naturally in 
the electroconvection of an annular suspended film that 
we have described. 

We next consider sheared, but non-rotating systems. 
RBC has been studied experimentally and theoreti- 
cally with an open through-flow. The through- 
flow is generally a weak Poiseuille flow with a very 
small Reynolds number. Its effects on RBC are well 
understood. In brief, the onset of convection is again 
suppressed, but the first instability is convective (i.e. 
it grows only downstream of a localized perturbation), 
rather than absolute. The resulting convection pattern 
drifts in the direction of the through flow. It is interesting 
that the clear distinction between convective and abso- 
lute instability is blurred in systems like ours, in which 
the "through" flow loops back on itself. As we have an- 
alyzed it, the onset of convection in our annular electro- 
convection is absolute. 

There have been some theoretical studies of 3D RBC 
in the presence of a plane Couette shear flow [|2| , another 
situation which is experimentally unrealizable. Here, lin- 
ear analysis reveals stability differences between trans- 
verse roll disturbances, those with axes are perpendicu- 
lar to the shear, and longitudinal roll disturbances which 
have axes parallel to the shear flow. Longitudinal-roll dis- 
turbances have identical stability properties to RBC in 



the absence of shear, and are always more unstable than 
the transverse- roll disturbances. Transverse-roll distur- 
bances, conversely, exhibit suppression, or added stabil- 
ity due to the shear. The onset Rayleigh number of trans- 
verse rolls is a monotonically increasing function of the 
shear Reynolds number, similar to what we find for 2D 
electroconvection. Furthermore, the critical wavenumber 
of the most unstable transverse disturbance was found 
to be a monotonically decreasing function of the shear 
Reynolds number, also as we observe. Transverse rolls 
(vortices, in fact) appear at onset in our 2D system, but 
would always be pre-empted by longitudinal rolls in a 3D 
RBC experiment, if it could be realized. 

IX. SUMMARY AND CONCLUSION 

In this paper, we have analyzed the linear stability of 
a weakly conducting, two-dimensional annular fluid to 
radially-driven electroconvection. We considered the ef- 
fects of an imposed Couette shear produced by the ro- 
tation of the inner edge of the annulus. We calculated 
the marginal stability boundary for the Rayleigh number 
TZ{a, V, TZe, m) as a function of the azimuthal mode num- 
ber m, for general radius ratio a, Prandtl number V and 
Reynolds number TZe. We also found the azimuthal trav- 
eling rate 7*/m of the marginally unstable modes, where 
7* is the imaginary part of the linear growth rate. We 
found a set of discrete values of a, V and TZe for which 
there are two adjacent azimuthal modes m which are si- 
multaneously unstable at onset. Two solutions schemes 
were employed: a local approximation which neglected 
some aspects of the coupling of fields and charges, and 
a more exact, fully nonlocal method. These methods 
agreed, except in some quantitative details. 

When there is no applied shear, so that TZe = 0, we 
found 7* = and the most unstable linear mode did 
not travel. As a — > 1, the marginal stability boundary 
for this case approached that for a laterally unbounded 
rectangular film, For non-zero shear, TZe ^ 0, the 
most unstable linear mode traveled in the direction of the 
rotation of the inner electrode and the onset of convection 
was suppressed, relative to the zero shear case. 

The shear also reduced the critical azimutal mode 
number tUc. Rigid rotation has no effect on the critical 
azimuthal mode number mc, again contrasting with the 
class of 3D rotating Rayleigh-Benard systems where 
increasing rotation rate favored higher critical azimutal 
mode number rric. 

We performed experiments on thin, suspended annu- 
lar films of smectic A liquid crystals, which are well de- 
scribed by our theoretical model. We measured the cur- 
rent through the film as a function of the driving voltage 
for a fixed rotation rate of the inner electrode. This data 
was used to determine the onset of convection and two 
of the material parameters needed to scale the data. We 
compared the relative suppression of the onset from the 
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experiments to that predicted by the Hnear theory and 
found good agreement. The spatial structure of the trav- 
ehng modes from the Unear theory resembled the weakly 
nonlinear vortex patterns observed experimentally, when 
viewed in a co-rotating reference frame. Q 

Due to the rich nonlinear aspects of this system, there 
are several interesting future directions. The system is 
sufficiently tractable that the weakly nonlinear regime 
should be accessible to theoretical investigation. It may 
also be interesting to simulate the nonlinear regime nu- 
merically. On the experimental side, we expect to be able 
to study the secondary bifurcations quantitatively, and 
make good connection to the nonlinear theory. It would 
be useful to determine the material parameters, partic- 
ularly the viscosity, by separate experiments. It should 
be possible to study regimes of higher TZe. Finally, it 
would be interesting to examine other flows, such as a 
time-dependent oscillatory shear, superposed on electro- 
convection. 
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APPENDIX A 



and {1, Bmn,Cmn, Dmn) IS the eigenvector correspond- 
ing to eigenvalue zero for each Pmn- We tabulated the 
Pmn at each radius ratio a using a multiple root finding 
routine. |30| The determinant was found using an SVD 
algorithm, and the eigenvectors by backsubstitution. |3^] 
It was necessary to rescale the third and fourth columns 
of M before taking the determinant to avoid underflows 
and overflows. 



thogonal set with orthogonality condition 



The functions defined by Eq. |A2 form a complete, or- 

(A4) 



The normalization Mcm -n is given by 
where 

It is convenient to define a function V„i-n by 
SO that C„i;„(r) = Um-n{r) + V„i;„(r) and 



D,D 



(A5) 
(A6) 
(A7) 



f3lni^m-n-Vrn,n)- (A8) 



1. Expansion functions for the stream function 

Since the stream function is constrained by rigid 
boundary conditions aX r = Vi and r = Tq and obeys 
Eq. 4.11. it can be expanded in the eigenfunctions of the 
square of the Laplacian operator. Hence, the eigen- 
functions sought are defined by the eigenvalue relation 



2. Expansion functions for the potential 

The potential perturbation satisfies homogeneous 
boundary conditions at the inner and outer edges of the 
film and obeys Eq. 4.12, so that we can use an expansion 



in terms of the eigenfunctions of the Laplacian operator. 



D^D 



(Al) 



mlWra;l i 



(A9) 



The boundary conditions are Cm 
and DCmni'Ti) = DCm-n{ro) = 0. 
of Eq. Al are 



The desired solutions 



(A2) 



where Jm and Ym are the familiar Bessel and Weber func- 
tions of order to, and Km are the modified Bessel 
functions of order to. The parameters Pmn are successive 
solutions of the secular equation det M(/3) = 0, where 



M(/3): 

JmiPn) YmiPn) imiPn) KmiPn) 

JmiPro) Ym{(3ro) ImiPro) KmiPro) 
Jm-liPn) Ym^liPri) Im^iiPn) -Km-liPn) 
Jm-iifiro) Ym-l{l3ro) Im-l{l3ro) -Km-liPVo) 



which has solutions 

r) . (AlO) 
solve the secular equation 



The parameters % 
det N(x) =0, where 



N(x) 



Jmixn) Ymixn) 

Jmixro) Ymixro) 



(All) 



and (l,6„m) is the eigenvector corresponding to eigen- 
value zero for each Xmn- These constants wer e fou nd by 
routines similar to the ones used in Section A 1 . The 



(A3) 



functions defined by Eq. AlO satisfy the orthogonality 
condition. 
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where the normahzation J^rp^.„ is given by 
1 



(A12) 
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